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Abstract 

Sensitivity to initial conditions in the coherent noise model of biological evolution, introduced 
by Newman, is studied by making use of damage spreading technique. A power-law behavior 
has been observed, the associated exponent a and the dynamical exponent z are calculated. 
Using these values a clear data collapse has been obtained. 
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It is clear that there is an increasing interest in extended dynamical systems exhibiting 
avalanches of activity, whose size distribution is scale-free. Some examples of such systems 
can be enumerated as earthquakes rice piles 0], extinction in biology Q], evolving 
complex networks ^, etc. However, there is no unique and unified theory for such systems. 
One of the possible candidates is the notion of self-organize criticality (SOC) l^. In SOC 
models, the whole system is under the influence of a small driving force that acts locally 
and these systems evolve towards a critical stationary state having no characteristic spatio- 
temporal scales. On the other hand, it is known from literature that another kind of 
simple and robust mechanism is available in order for producing scale-free behavior. This 
mechanism is based on the notion of external stress coherently imposed on all agents of 
the system. Since the model does not contain any direct interaction among agents, it does 
not exhibit criticality. Nevertheless, it yields a power-law distribution of event size. These 
so-called coherent noise models have been firstly introduced to describe large-scale events 
in evolution, but then they were used as a model of earthquakes I? I and its aftershock 

ft n 

properties [q] as well as aging phenomenon in the model ^9i] have been analysed. 

The coherent noise models can be defined as follows: Firstly we consider a system, 
which has agents. Each agent i has a threshold Xj against external stress 77. The 
threshold levels are chosen randomly from some probability distribution pthresh{x)- The 
external stress is also drawn randomly from another distribution Pstressiv)- -^"^ agent 
is eliminated if it is subjected to the stress r] exceeding the threshold for this agent. 
Algorithmically, dynamics of the model can be given by three steps: (i) at each time step, 
generate a random stress r] from Pstressiv)^ eliminate all agents with Xi < i] and replace 
them by new agents with new thresholds taken from pthresh{x), (ii) select a small fraction 
/ of the agents at random and assign them new thresholds, (iii) go back to (i) for the 
next time step. It is worth mentioning that step (ii) corresponds to the probability for 
the / fraction of the whole agents of undergoing spontaneous transition, which is a step 
necessary for preventing the model from grinding to a halt . 

In the present work, we focus on the sensitivity to the initial conditions of the coherent 
noise models using damage spreading technique. This technique can be thought as a 
method which is borrowed from dynamical systems theory in the following sense: If we 
consider two copies of the same one-dimensional dynamical system starting from slightly 
different initial conditions and follow their time evolution, we can define the sensitivity 
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function 



^(t)^ lim = e'* (1) 

to quantify the effect of initial conditions. Here, Ax(0) and Ax{t) are tlie distances 
between two copies at t = and t respectively, and A is the Lyapunov exponent. If A > 
(A < 0), the system is said to be strongly sensitive {strongly insensitive) to the initial 
conditions. For the marginal case, where A = 0, the form of the sensitivity function could 
be a whole class of functions. For the low- dimensional discrete dynamical systems, this 
form is found to be a power-law 

m-t- ■ (2) 

The a > and a < cases correspond to weakly sensitive and weakly insensitive to the 
initial conditions Q| . For the high- dimensional dynamical systems, like the Bak-Sneppen 
model or like the one that we discuss in this work, the same analysis could be performed 
using the Hamming distance instead of the sensitivity function. As in the case of low- 
dimensional systems, one can classify the sensitivity by looking at the behavior of the 
Hamming distance 

D{t) = l^T.\^^'\t)-xf{t)\ , (3) 
1=1 

where x\^^ and xf'^ are two slightly different copies of the system under consideration. In 
this way, up to now, various variants of the Bak-Sneppen model have been analysed and 
related exponents are calculated jlll . 12, ^ Q, Q|. Our aim here is to make the same 



analysis for the coherent noise models to numerically obtain the dynamical exponents of 
the model. 

To proceed further let us define the procedure. We start the simulations from a uniform 
threshold distribution, which means that the system starts from an initial event of infinite 
size that spans the whole system. This is considered to be the first copy of the system 
(namely, x[^^). In all simulations, we use the exponential distribution for the external 
stress 



Pstressiv) = a ^exp ( -^j (a > 0) , (4) 
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and the uniform distribution pthreshi^) (0 < x < 1) for the threshold level. The second 
copy (namely, x\^^) is generated by exchanging the values of two randomly chosen sites. 
Then the dynamics of the model, as explained above, is applied for both copies of the 
system using always the same set of random numbers to update them. To investigate 
the properties of sensitivity to initial conditions of the model, we follow the temporal 
evolution of the Hamming distance given in Eq.(3). In all cases, we use a large number 
of realizations to reduce the fluctuations and the given results, namely {D{t)), are the 
ensemble averages over these realizations. Prom Fig. 1 it is easily seen that the damage 
spreads as a power-law, indicating a weak sensitivity to initial conditions. Obviously, 
the power law growth {D{t)) ~ with a = 0.95 ± 0.04 is followed by a plateau with a 
constant value, starting at a certain time depending on the system size. This saturation is 
due to the fact that by this measure one cannot identify that both copies have converged 
to the same random sequence. 

The second important exponent, called as the dynamical exponent z, comes from the 
scaling of t{N), which is defined to be the value of t where the power-law increasing 
part crosses over onto the saturation regime. More precisely, r is the value of t at the 
intersection point of two straight lines, one of which comes from the power-law curve and 
the other from the constant plateau. Then, the scaling t{N) ~ defines the dynamical 
exponent z, which is obtained from Fig. 2 as z — 0.96 ± 0.02. 

As a final step, we analyse the finite size scaling behavior of the model. To accomphsh 
this task, we define the normalized Hamming distance as 

and numerically verify that it obeys the scaling ansatz 

D(iV,i) = Ar^F(^) , (6) 

with /3 = 0.91, which comes from P — ja and ^ — z. A clear data collapse can easily 
be seen in Fig. 3. It is worth noting that the results presented here do not depend on 
the choice of the numerical value of the model parameter a and/or the fraction /. In our 
simulations we use a = 0.001 and different / values for each system size in order to assure 
only one agent is eliminated in step (ii) to reduce the computational time. 
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FIG. 1: Time evolution of the normalized Hamming distance for five different system sizes. For 
the largest one (A^ = 32000) the exponent of the power-law growth is estimated as a = 0.95zb0.04 
The number of realizations used for averaging is 200 for N = 32000 and 400 otherwise. For a 
better visualization logaritmic binning is employed to all data sets. 

Summing up, we have studied the sensitivity to initial condition properties of the 
coherent noise models using damage spreading technique. We found that the model 
exhibits weak sensitivity to initial conditions, a property which is common for other high- 
dimensional dynamical systems like Bak-Sneppen model and its variants as well as one- 

and two-dimensional discrete systems like logistic map families. The numerically obtained 
values of the power-law exponent a and the dynamical exponent z are appeared to be 
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FIG. 2: Double log plot of r versus N for each curve given in Fig. 1. From this scaling we 
estimate the dynamical exponent as z = 0.96 it 0.02. 

different from those of the Bak-Sneppen model, which signals out the discrepancy between 
the universality classes of these models. 
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FIG. 3: Data collapse of the finite size scaling given in Eq.(6) for the same system sizes used 
in Fig. 1 and Fig. 2. 
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